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Networks 
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Abstract — This paper considers the problem of detecting topol¬ 
ogy variations in dynamical networks. We consider a network 
whose behavior can be represented via a linear dynamical system. 
The problem of interest is then that of finding conditions under 
which it is possible to detect node or link disconnections from 
prior knowledge of the nominal network behavior and on-line 
measurements. The considered approach makes use of analysis 
tools from switching systems theory. A number of results are 
presented along with examples. 

I. Introduction 

Recent years have witnessed a growing interest towards 
networks of dynamical systems. There is in fact a trend to 
build modern infrastructures as large-scale networks, which 
are possibly geographically distributed [l]-[4]. Networks of 
dynamical systems also arise and play a fundamental role 
in vehicle formation, cooperative robotics, surveillance and 
environment monitoring, to name a few [5]-[8]. 

In many situations, the network behavior is determined or 
strictly related to its underlying topology. This is the case, 
for instance, in consensus, coordination and synchronization 
problems, where the dynamical system that describes the 
evolution of the network is related to the structure of the 
graph that models the interaction among the various networks 
components [9], [10]. On the other hand, also the problem 
of inferring the network topology from observations of the 
network behavior is of paramount importance, and it is the 
objective of this paper to explore such a topic. 

Variations of the network topology can have a major impact 
on stability and/or performance. For example, in consensus¬ 
like networks a link disconnection may slow down conver¬ 
gence or even destroy agreement when the graph connectivity 
is lost [7], [9]. More importantly, variations of the topology 
may affect the network secure and reliable operation. In fact, 
the strong interdependency among the various elements of the 
network is such that a failure in one part of the communication 
infrastructure can rapidly create global cascading effects. This 
issue is amplified by the fact that failures in the communication 
infrastructure can be caused not only by equipment failures or 
human errors but also by intentional attacks [11]-[13] 

This paper considers the problem of inferring variations 
of the network topology from observations of the network 
behavior. Specifically, we consider a network whose behavior 
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can be represented via a linear dynamical system. The problem 
of interest is then that of finding conditions under which it 
is possible to detect a node or link disconnection from prior 
knowledge of the nominal network topology and measure¬ 
ments of the network state or a subset of it. Contributions 
to this topic have been recently proposed. In [14], [15], the 
authors address the problem of detecting single and multiple 
link failures in a multi-agent system under the agreement 
protocol. A notion of distinguishable flow graphs is introduced 
and sufficient conditions for achieving distinguishable dynam¬ 
ics are stated in terms of inter-nodal distances. In [16], the 
authors investigate the problem of detection and isolation of 
link failures by exploiting the presence of discontinuities in 
the derivatives of the output responses of a subset of nodes. 
It is worth noting that the problem of inferring variations 
of the network topology can also be addresses by means of 
topology identification algorithms [17]-[19]. However, iden¬ 
tification algorithms do not assume prior knowledge of the 
nominal network topology, which is possible in many practical 
circumstances, and, as such, do not take full advantage of such 
extra information, which may be crucial for achieving early 
detection of stability and/or performance losses. 

The approach taken in this paper makes use of analysis tools 
from switching systems theory. Specifically, networks with 
switching topology can be naturally modeled as a switching 
system, where the switching signal determines the current 
network configuration (operating mode). Thus, the problem 
of detecting a node/link disconnection can be naturally cast 
as the problem of determining under what conditions the 
operating mode of the system can be uniquely reconstructed 
from observations. In the relevant literature, this problem 
is known as the discemibility, distinguishability or mode- 
observability problem [20]-[26]. 

For linear systems, discemibility can be fully characterized 
through simple algebraic conditions. In fact, it is completely 
characterized by the eigenspace components related to the 
various operating modes of the system. These conditions are 
generally difficult to refine because the dynamics related to 
the various operating modes of the system need not be related 
with one another. In the present case, however, the situation 
is different because the dynamics resulting from a node or 
link disconnection can be related with the nominal one via 
interlacing theorems [27]. Moreover, for several graphs of 
practical relevance, such as complete, ring, path and grid 
graphs, an explicit expression for the eigenspace components 
is available. 

By exploiting these features, we provide necessary and 


sufficient conditions for detecting topology variations for both 
the cases of node and link disconnections. These conditions 
are based on simple algebraic tests, which can be easily 
checked numerically as well as analytically whenever an 
explicit expression for the eigenspace components turns out 
to be available. While the analysis is mainly oriented towards 
a theoretical characterization of the detection problem, the 
results also provide several insights on how detection can be 
addressed in practice, as well as guidelines for the develop¬ 
ment of sensor placement algorithms. 

The remainder of this paper is as follows. In Section II, we 
describe the framework of interest and formulate the detection 
problem. In Section III, the main results of the paper are given. 
Connections with least-square identification are established in 
Section IV. Finally, Section V ends the paper with concluding 
remarks. For convenience, the proofs are reported in the 
Appendix. 

II. Framework and Problem Formulation 

We consider a network of n nodes, whose topological 
structure is represented by an undirected graph Q := {y,£), 
where V := {1,2,... ,n} denotes the node set and £ CV xV 
denotes the edge set. We assume that the network behavior 
can be represented via a linear dynamical system 

X = ( 1 ) 

where x G M" denotes the network state; Xi G K”, i G V, 
denotes the state of the i-th network node; <I> G is the 

matrix that determines the network behavior. We assume that 
$ = $' and that (j)ij 7^ 0 if and only if {i,j) G £, where (j)ij 
denotes the (i,j)-th entry of <I>. 

As a relevant example, consider a classical agreement 
problem in a network of continuous-time integrators with 
local dynamics Xi = Ui, which implement a linear consensus 
protocol with unitary weights, 

^ {Xj - Xi) ( 2 ) 

jeAfi 

where Afi denotes the set of neighbors of node i. This gives 
rise to the linear system x = —Lx, where L denotes the graph 
Laplacian induced by Q. The system is therefore in the same 
form as ([T]i with $ = —L, i.e., with (j)ij = 1 for j ^ i and 
(l)ij = -|M| for j = i. 

Remark 1: Although this paper is only concerned with 
networks whose topological structure is represented by an 
undirected graph, most of the conclusions can be extended 
to directed graphs as well. ■ 

A. Problem formulation 

We regard the pair (Q, <I>) as representative of the nominal 
behavior of the network. The problem of interest is then that 
of hnding conditions under which it is possible to detect a 
variation of the network topology via: i) knowledge of the 
nominal matrix $; and ii) measurements of 

(3) 


where S K" denotes the i-th versor, and AA C V denotes 
the set of nodes whose state is available for measurements. 

A variation in the network topology is specified by means of 
a pair {Q, $), where Q describes the novel topological structure 
and $ describes the novel network behavior, i.e., 

X = ^x ( 4 ) 

In particular, we assume that Q is an undirected graph defined 
as Q := {V,£), where f C f; $ = $' with fij 7^ 0 if and 
only if {i,j) G £, where fij denotes the (i,j)-th entry of 
$. As detailed hereafter, {Q, $) captures several scenarios of 
practical relevance. 

1) Link disconnection without dynamics reconfiguration: 
suppose that the link disconnection affects the nodes i,j G V. 
Then Q is characterized hy £ = £ \ {{i,j), {j, i)}, and 

|> = $ - (e^e' -f Cjc'^ (5) 

In words, the network dynamics remains unchanged with the 
exception that fiij = fiji = 0. 

2) Link disconnection with dynamics reconfiguration: this 
scenario is the same as the previous one with the exception 
that, in addition to having 0^ = fiji = 0, a variation occurs 
also in fin and fijj. For example, in the linear consensus 
problem described above, one has fin = fin + l and fijj + 1, 
and 

$ = $ -p ae'i A- tje'j - eie'j - Cjc'i (6) 

3) Node disconnection without dynamics reconfiguration: 
Suppose that the node disconnection affects the node i G V. 
Then Q is characterized hy £ = £ \ {{i,j), (j, *); j G Afi}, 
while 

i = 4* - ^ {cie'j + eje[) (7) 

j&Ni 

4) Link disconnection with dynamics reconfiguration: this 
scenario is the same as the previous one with the exception 
that a variation occurs also in fin. In the linear consensus 
problem described above, one has fin = 0, and 

$ = $ -f ^ (cic' -f Cje'j - ae'j - Cje'i) (8) 

j&Ni 

III. Main Results 

In this section, we first introduce a notion of discernible 
networks. We then present the main results of this paper 
and establish a number of connections with several graphs 
of practical interest. To begin with, notice that the problem 
of detecting a node or link disconnection from measurements 
can be casted as the problem of finding conditions under 
which $ and $ do not give rise to the same dynamics. 
Networks satisfying this property can be therefore referred to 
as discernible. 

We formalize these concepts. 

Definition 1: A dynamical network is said to described by 
the pair {Q, $) if Q is the graph describing the network 


y = Mx, M = col{ei, i G AA} 


topology and the network behavior obeys x = where x is 
the network state. ■ 

Definition 2: Consider two dynamical networks described 
by {Q, $) and {Q, $), respectively. The networks are said to 
be indiscernible with respect to the state a;o if = e**a;o 

for all t £ ]R>o Otherwise, they are said to be discernible. We 
denote by I the set of states for which {Q, $) and {Q, l>) are 
indiscernible. ■ 

Definition 3: Given a matrix M as in Q, the networks are 
said to be M-indiscernible with respect to the pair of states 
(a:o,a;o) if Me'^*Xo = Me^*xo for all t £ M>o- Otherwise, 
they are said to be M-discernible. We denote by X{M) the 
set of pairs of states for which [Q, $) and (^, $) are M- 
indiscernible. ■ 

Both discernibility and M-discernibility can be viewed as 
particular observability problems, which can be addressed by 
looking at the parallel interconnection of i and x = $a;. 

For instance, discernibility is equivalent to the observability of 
the pair (A, F), where A = diag{$, $} and F = [J — /], over 
the set X = {{w, : w = ^}. On the other hand, M- 

discernibility is equivalent to the (standard) observability of 
(A, F), where A = diag{$, !>} and F = [/ — /]. The latter is 
the classical condition for reconstructing the active mode of a 
switching linear system from output measurements [20], [21]. 
One sees that both discernibility and M-discernibility depend 
entirely on the eigenspaces of $ and $, which are in general 
difficult to analyze. With respect to the case of switching 
systems, however, the analysis here considerably simplifies 
since $ and $ comes with a symmetric structure. For clarify, 
we address the cases of discernibility and M-discernibility 
separately. 

A. Discernibility 

We first consider the discernibility problem. Notice that 
since $ and $ are symmetric, there exist orthonormal matrices 
S and S such that 

$ = SAS', !> = SKS' (9) 

with A and A diagonal matrices. Let spec($) denote the 
spectrum of <i>. Moreover, for any A £ spec($), let /i(A) 
and V (A) denote its multiplicity and eigenspace, respectively. 
Finally, let S'(A) be the set of columns of S that generate 
L(A), i.e., V(A) = span(S'(A)), where span denotes the linear 
span. We then have 


^ e"‘5(A)5'(A) 

(10) 

AGspec($) 


^ e"‘5(A)5'(A) 

(11) 

AGspec($) 



From the above expressions, it is straightforward to draw the 
following conclusions: 

i) If spec(<i)) n spec(T>) = 0, then X = {0}, i.e. the only 
indiscernible state is the zero state. This is obviously the 
smallest indiscernibility set that one may have. 


ii) Nonzero indiscernible states exist if and only if there 
exists some eigenvalue A common to <!> and such 
that V{X) n L(A) 7 ^ 0, or, equivalently, such that 
rank{[S'(A) 5(A)]} < p.(A) -I- fi{X). 

Let 'I'(A) be a matrix whose columns form an orthonormal 
basis of L(A) n L(A). Hence, the set X of states for which 
{0, $) and {Q, T>) are indiscernible is given by 

X = \^x : X € span ('I'(A), A £ spec($) n spec(l>))} (12) 

In view of the above considerations, it is interesting to 
investigate under which circumstances the variation in the 
network topology may lead to a new system matrix $ sharing 
eigenvalue/eigenvector pairs (A, x) with the original system 
matrix $. Clearly, this amounts to searching for necessary and 
sufficient conditions for the existence of pairs (A, x) such that 

= Xx 

With respect to the four scenarios of interest described in 
Section II-A, the following results can be stated, which show 
that indiscernible states can be readily inferred by inspection 
of the components of the eigenvalues of $. 

Theorem 1: (Link disconnection without dynamics recon¬ 
figuration). Consider a disconnection of link {i,j) as in (|5]l. 
Then, the networks are indiscernible with respect to a state 
X £ L(A), with A £ spec(‘l>), if and only if Xi = xj =0. ■ 

Theorem 2: (Link disconnection with dynamics reconfig¬ 
uration). Consider a disconnection of link {i,j) as in (|6]l. 
Then, the networks are indiscernible with respect to a state 
X £ L(A), with A £ spec(‘l>), if and only if Xi = xj. m 

Theorem 3: (Node disconnection without dynamics recon¬ 
figuration). Consider a disconnection of node i as in ©. 
Then, the networks are indiscernible with respect to a state 
X £ L(A), with A £ spec($), if and only if 

Xi ^ ] 4^ij — ^7 ^ (13) 

ieW jeW 

If, in addition, we assume that fiij > 0 for any i j, then 
condition ([T3I) becomes Xi = 0 and Xj = 0 for any j € fifi. 

■ 

Theorem 4: (Node disconnection with dynamics reconfigu¬ 
ration). Consider a disconnection of node i as in Then, the 
networks are indiscernible with respect to a state x £ L(A), 
with A £ spec($), if and only if 

Xj = Xi, WjGjVi. (14) 

If, in addition, we assume that $ = —L (as in the linear 
consensus protocol (|2]i) and we consider a non-null Laplacian 
eigenvalue A, then condition ([T3l) becomes Xi = 0 and xj = 0 
for any j G Aft. m 

From the above results, it can be seen that, in general, 
node disconnections are easier to detect than link disconnec¬ 
tions. Similar considerations can be made to conclude that a 
link/node disconnection without dynamics reconfiguration is 
easier to detect as compared to a disconnection with dynamic 
configuration. 


As an example, consider again the case <i) = —L, where L is 
the graph Laplacian of Q. As well known, the all-ones vector 
1 is always an eigenvector of L associated to the eigenvalue 0. 
After a link/node disconnection with dynamic reconfiguration, 
the novel dynamic matrix $ will coincide with —L, where L 
is the Laplacian of the graph Q. Hence, the all ones vector 1 
will be an eigenvector associated to the 0 eigenvalue also for 
$ = —L. Hence, in this case, the stationary state a; = 1 turns 
out to be indiscernible for any link/node disconnection. This is 
consistent with the fact that a: = 1 satisfies the conditions of 
Theorem 2 and 4 for any i and any j, being all its components 
identical. On the contrary, when a link/node disconnection 
without dynamic reconfiguration occurs, a; = 1 is no longer an 
eigenvector of $ and, hence, is discernible (indeed x = 1 does 
not satisfy the conditions of Theorem 1 and 3). If we further 
assume that the graph Q is connected, then 1 turns out to be 
the unique eigenvector with eigenvalue 0. Then, the second 
part of Theorem 4 can be exploited to conclude that, with 
the exception of x = 1, the two cases of node disconnection 
with or without dynamic reconfiguration give rise to the same 
indiscernible eigenvectors. 

B. M-Discernibility 

We now turn the attention to the problem of detecting a 
topology variation from observations of an output vector y = 
Mx. In this respect, notice preliminarily that a state x G I 
for which {Q, $) and {Q, $) are indiscernible does always 
generate indiscernible output trajectories. Then, if we define 
the set Xp = {(x,x) : x G I}, we have that Ip C I(M) 
irrespective of the choice of the output matrix M. On the other 
hand, by letting M = I (i.e., by observing all the network 
nodes) we clearly have Ip = I(/), where / stands for the 
identity matrix of an appropriate dimension. Hence, a first 
important problem is how to choose the matrix M so that Ip = 
I{M). This amounts to asking where sensors nodes should 
be placed in order to guarantee that discernibility implies M- 
discernibility. When such a condition holds, we say that the 
matrix M ensures output discernibility. Since we have 

e^^MS{\)S'{X) 

AGspec(^) 

Me** = Y e^*MS{X)S'{X) 

AGspec(^) 

the following result follows at once. 

Theorem 5: Consider an observation vector y as in 0. 
Then, condition Ip = I(M) holds if and only if the following 
conditions are satisfied; 

(i) rank{MS'(A)} = /r(A) for any A G spec($) \spec(i>); 

(ii) rank{M;S(A)} = fl{X) for any A G spec($) \spec($); 

(iii) rank {M [S'(A) >S(A)|} = rank { [S{X) 5(A)] } for any 
A G spec($) n spec($). 

■ 

Notice that condition (i) amounts to requiring that all the 
states belonging to V (A), with A eigenvalue of $ but not of $, 
are observable from the output y = Mx. The same property is 


required by condition (ii) for all the eigenvalues of $, which 
are not eigenvalues of $. Finally, condition (iii) amounts to 
requiring that, for any eigenvalue A that is shared by $ and 
$, and for any (x,x) G V{X) x C(A), one has Mx = Mx if 
and only if x = x. 

Building upon Theorem 5, an expression for I(M) can be 
given. To this end, for any A G spec(<I>)\spec($), let K{X, M) 
be a matrix whose columns form a basis of the linear space 
{(x,0) G : X G y(A) and Mx = 0}. Let K{X,M) 
be defined in a similar way with respect to A G spec(<I>) \ 
spec($). Finally, for any A G spec($) nspec($), let T(A, M) 
be a matrix whose columns form a basis of the linear space 
{(x,x) G : X G V{X),x G C(A), and Mx = Mx}. 
Then, we have 

I(M) = span {T(A, M),X G spec($) n spec)®)} 

Uspan{Ar(A,M), A G A G spec($) \spec(i>)} 
Uspan{Ar(A,M), A G A G spec($) \spec($)} 

Theorem 5 provides interesting insights on the number of 
sensors needed so as to have output discernibility. Consider, 
for example, the ideal situation in which $ and $ are dis¬ 
cernible from all the states, i.e., I = {0}. Then, condition 
(iii) becomes 

rank {M [S{X) 5(A)] } = p.{X) + fl{X) 

for any A G spec(<I>)nspec($). Notice also that the rank in the 
left-hand side cannot exceed rank{M}, which, in turn, is equal 
to the number of measured nodes. Then, we can conclude that, 
in order to have output discernibility, one needs a number 
of sensors at least equal to the maximum among y,{X) for 
A G spec($) \spec($), p,{X) for A G spec($) \spec($), and 
fj,{X) + fix) for A G spec($) n spec(#). Note that this is 
just a lower bound, since Theorem 5 does not exclude that a 
larger number of sensors may be needed. Nevertheless, such 
considerations indicate that, similar to what happens when 
standard observability is addressed [28], the number of nodes 
that should be available for measurements increases with the 
multiplicity of the eigenvalues. 

As for the sensor placement, one can see that in order to 
satisfy condition (iii) the sensors i € A4 must be positioned so 
that the rows of the matrix [S'(A) 5(A)] corresponding to the 
indices i € A4 contain at least one non-zero minor of order 
rank{[5'(A) 5(A)] }. Analogous considerations can be given 
for conditions (i) and (ii). 

In particular, condition (ii) becomes tricky when, starting 
from a connected graph Q, a topology variation gives rise 
to multiple connected components in the graph Q (notice 
that this always happens in the case of node disconnection). 
Specifically, let Q consist of N mutually disjoints components 
0^,..., , and let be the set of nodes belonging to 

(clearly l-^^l = Then, as well-known, 

N 

spec(^) = (J spec(^^) 


and, in addition, for any A G spec(^^) there exist eigenvectors 
X G t^(A) such that a;^ 7 ^ 0 if and only if i S A/”^. As a 
consequence, it is immediate to verify that condition (ii) can 
be satisfied only by placing sensors in each one of the mutually 
disjoints components Q^,..., . 

Clearly, this latter requirement can be quite restrictive in 
practice. For instance, this implies that one can have output 
discernibility with respect to any possible node disconnection 
only by placing a sensor in each network node. Hence, instead 
of requiring complete output discernibility, in many situations 
it may be of interest to restrict the attention only to some 
of the connected components of the graph Q. This can be 
done in a straightforward way by considering in condition 
(ii) only the eigenvalues and eigenvectors pertaining to the 
connected components of interest. For example, in the case of 
disconnection of node i, one can restrict the attention to the 
component with node set A/^\ {i} by excluding from condition 
(ii) the eigenvector Ci pertaining to the trivial component {i}. 

IV. A LEAST-SQUARES CRITERION EOR DETECTION OE 
TOPOLOGY VARIATIONS 

In the previous section, we have provided conditions under 
which it is theoretically possible to detect a topology variation 
by observing the evolution of the state Xi in a subset AT of 
the network nodes JV. From a practical point of view, this can 
be done by resorting to a least-squares criterion, as detailed 
hereafter. 

Suppose that, starting from time tg, a certain number, say 
N, of samples of the output vector y are collected at the 
time instants t = tg + kT for k = 0,..., A^ — 1, where 
T G K.>o- In particular, to account for the possible presence 
of a measurement noise, let each sample Zk be of the form 

Zk = y{tQ + kT) + Vk ( 15 ) 

with Vk an unknown but bounded discrete-time noise signal. 
We assume that an upper bound Ey on the energy of the 
sequence {vk} is known, i.e., where 

II • II stands for Euclidean norm. Hereafter, the vector of all the 
collected samples will be denoted by 

Zpf = co\{zo,...,zn-i) 

Remark 2: It is worth noting that (fTSl l amounts to making 
use of synchronous measurements. While this hypothesis may 
be restrictive in some cases, there are many applications where 
the measurement devices are equipped with global positioning 
system (GPS) units. This is the case, for instance, in many 
smart grid applications where Phasor Measurement Units are 
sampled from widely dispersed locations and synchronized via 
a common GPS reference [2]. ■ 

Let On and On denote the sampled-data observability 
matrices associated with Q and Q, respectively. Clearly, we 


have 



( ^ \ 


( ^ \ 





On ■= 


, On ■■= 







( 16 ) 


Notice now that, when the state evolution is generated by the 
nominal network ((/, $), the sampled outputs are of the form 
Zn = ONx{to) + Vn where Vn = col(t;o,..., uat-i) and 
x{to) is the (unknown) state at time to. Then, the least-squares 
cost function 

iriZN) = min \\Zn - On x\\ 

provides a quantitative measure of how close the observed 
output behavior is to the nominal ones. In fact, whenever the 
output samples Zn arises from the nominal network, we have 
7r(ZAr) < Ey. 

Similarly, any output behavior generated by {Q, $) leads to 
sampled outputs of the form Zn = ONx(tQ) + VN and, hence, 
the least-squares cost function 

Tr{ZN) = min \\Zn - On x\\ 

provides a quantitative measure of the distance between the 
observed outputs and the set of behaviors associated with the 
modified topology. 

Then, by computing the quantities 7r{ZN) and tt^Zn), the 
following conclusions can be readily drawn; 

(a) when ^{Zn) > Ey, the output samples are not consistent 
with the nominal network; hence we can conclude that a 
variation from the nominal behavior has occurred. 

(b) when ^{Zn) < Ey and ^{Zn) > Ey, the output samples 
are consistent only with nominal behavior; hence we can 
exclude the variation associated with ((/,<&). 

(c) when both ^{Zn) < Ey and ^{Zn) < Ey, we cannot 
conclude since the sampled outputs are consistent with 
both the nominal and the modified behavior. 

Of course, the considered framework can be easily extended so 
as to account for different possible topological variations. In 
particular, the detection of a topological variation as in case (a) 
requires only the computation of the cost function associated 
with the nominal behavior. On the other hand, a “validation” 
test as in case (b) requires, in general, the computation of one 
cost function for each possible variation. Further, also in case 
(a), computation of the other cost functions can be useful in 
order to possibly identify the topological variation. 

As for case (c), this corresponds to the situation in which 
the information contained in sampled outputs is not sufficient 
to conclude on the underlying topology. Clearly, this case 
may arise when the measurement noise is sufficiently large 
so as to mask the data information. However, case (c) is also 
inherently connected to the concept of M-indiscernible states 
as previously defined. In fact, let I{M) and I{M) be the 
projections on the first and, respectively, last n components of 






the set I{M), i.e., 

T{M) = {cc G K" : e M” such that (a;, x) G I{M)] 

X{M) = {x G K" : 3a: G M" such that (a:, x) G X{M)} 

Then, it is an easy matter to see that when the output behavior 
is generated by the modified network starting from state 
x(<o) G X{M), one has tt{Zn) < Ey and 7 r(Zjv) < Ey since 
both TT{ONx{to)) = 0 and TT{ONx{to)) = 0. In this respect, 
it is worth noting that although X{M) is defined with respect 
to an ideal situation (i.e., assuming to measure the noise-free 
continuous-time evolution of y = Mx), it turns out that such 
a set plays a fundamental role also in the practical situation 
of sampled outputs affected by measurement noises, provided 
that the sampling is non-pathological [29]. 

Lemma 1: Let the state trajectory be generated by the 
modified network (^, $). Furthermore, suppose that for any 
A, A G spec(<l>) Uspec($) with A 7^ A the following condition 
holds 

Im(A — A) 7 ^ for G Z \ {0} 

whenever Re(A — A) = 0 . (17) 

Finally, let N > 2n. Then, case (c) can occur only if 

d{x{to),X{M)) <aEy (18) 

where a is a suitable positive constant and stands for 

point-set distance. ■ 

In view of Lemma 1, it can be seen that, when the state 
a;(fo) is far enough from the set of M-indiscernible states (in 
the sense that condition (ITsT l does not hold) and the sampling is 
non-pathological, then the sampled outputs provide sufficient 
information for detecting that a topological variation has 
occurred. Bounds on the constant a can be found, as in [25], in 
terms of the cosine of the smallest non-null angle between the 
linear subspaces spanjOAr) and span((!)Ar). Notice finally that 
condition (ITtI i is nothing but the well-known Kalman-Bertram 
criterion for the observability of sampled-data systems applied 
to the pair (F, A). 

Notice that a result analogous to Lemma 1 can be derived 
also in the case of output behaviors generated by the nominal 
network (0, $), with the set X{M) replaced by X{M). 

V. An Example 

In order to illustrate some of the previous results in an easy 
manner, we consider the simple, yet non trivial, case of a 
linear consensus protocol with unitary weights over the n- 
dimensional path graph G = Pn with dynamics reconfigura¬ 
tion, which is standard in consensus-like algorithms. We focus 
on the discernibility problem. Recall that the Laplacian L of 
Pn has eigenvalues A^^l = 2 — 2cos(TTk/n) and eigenvectors 

x^^^ = cosiy-Kki/n — 7rfc/2n) 

k G {0, l,...,n — 1},J G {1,2,...,n} 

where A^^l denotes the fc-th eigenvalue and x^^^ denotes the i-th 
component of the k-th eigenvector. This describes the nominal 
network {G, —L). 



Fig. 1. Indiscernible points for the path graph in case of link disconnections. 
Letting A := ■Kkiln — i^kl'ln and B := •Kkiln+i^kl'ln, one has cos(A) = 
cos(B) whenever A = ±m7r — A/2 and B = A + A, where m G N and 
A g [0, 7r). 


We consider first the case of disconnection of a link (i,j). 
Notice that we can restrict the attention to the situation where 
j = i + 1 since the case j = i — 1 is specular. In accordance 
with Theorem 2, the link disconnection is not detectable if and 
only if 

cos{'Kki / n — Trk/2n) = cosiirkj jn — 7rfc/2n) 

= cosiTikijn 3- 'Kkl2n) 

for some k G {0,1,..., n — 1} and i G {1, 2,... ,n — 1}. 
For convenience, let us define A := irki/n — Trkj2n and 
B := irki/n + 7rkj2n. By looking at the cosine function 
(cf. Figure one sees that discernibility is violated when 
A and B take the form A = ztmTr — A/2 and B — A + A, 
where to G N := {0,1,2,...} and A G [ 0 , 7 r). The reason 
for constraining A to be less than tt comes from the fact that 
B = A + 7rk/n < A + tt since k < n. This in particular 
excludes indiscernible points located at the zeros of the cosine 
function as well as indiscernible points of periodicity of 27r 
or higher. Combining the previous expressions, we then have 
A = Trk/n Thus discernibility is violated whenever exist 
k,i,n and to such that A = 7 rfci/n — 7 rfc/ 2 n = Erfi'K — 'Kkl2n 
or, equivalently, 

ki = nm (19) 

where n > 1. The following conclusions can be drawn: i) A 
trivial solution to is given by fc = to = 0. As previously 
noted, this corresponds to the fact that the consensus state x = 
1 turns out to be indiscernible for any link/node disconnection; 
ii) Nontrivial solutions to ( [T9] ) also exist. Simple examples are 
{k,i,n,m) = (2,4,8,1) and {k,i,n,m) = (8,5,10,4); iii) 
Apart from the case where k = m = 0, there is no solution 
to dT^ when i G {l,n — 1}. In fact, if j = we obtain k = 
nm. This has no solutions when k,m > 0 since fc < n. If 
instead / = n — 1 we obtain fc(n — 1) = nm, or, equivalently, 
k — m = k/n. This has no solutions when k,m > 0 since 
fc — TO is integer, whereas k/n cannot be integer since k < n. 







Point iii) is interesting since it indicates that, apart from the 
stationary state a; = 1, a link disconnection can always be 
detected if it involves one of the endpoints of the graph. The 
latter situation can also be viewed as a node disconnection, 
and, in fact, it is straightforward to verify that, apart from the 
stationary case, node disconnections are always detectable. To 
see this, notice that, in accordance with Theorem 4, a node 
disconnection is not detectable if and only if there exist an 
index i such that 

cosii^ki! n — 7rfc/2n) = cosij^ki j n + 'nkj^ln) 

= cosiyirki / n — 37rfc/2n) 

By the point iii) above, one can restrict the attention to the 
case where i G {3,4,..., n — 2} since the cases i G {1, 2} 
and i € {n — l,n} do involve the endpoints of the graph. Let 
A and B be as before, and let C := irki/n — 37rfc/2n. Thus, 
there must exist three points, namely A, B and C where the 
cosine function takes on the same value. Apart from the case 
k = 0, this is not possible since B = C + 2TTk/n < C + 2 tt. 
Hence, we conclude that, apart from the stationary case, a 
node disconnection is always detectable. 

It is worth mentioning that a very similar analysis could be 
carried out with respect to the grid graph since its eigenspace 
is completely determined by the eigenspace of the path graph 
[28]. 

VI. Conclusions 

In this paper, we have addressed the problem of detecting 
topology variations in dynamical networks, considering both 
the cases of node and link disconnections. The results show 
that the detection problem can be characterized through simple 
algebraic conditions, which depend on the eigenspace com¬ 
ponents related to the nominal and faulty operating mode of 
the network. While the analysis is mainly oriented towards 
a theoretical characterization of the detection problem, the 
results also provide several insights on how detection can be 
addressed in practice, as well as guidelines for the develop¬ 
ment of sensor placement algorithms. 

Appendix 

Proof of Theorems 1-4. Let x be an eigenvector of $ with 
eigenvalue A, i.e., x is such that = \x. Then, Theorem 1 
readily follows from the fact that, when a link disconnection 
without dynamics reconfiguration occurs, we have 

^x = ^x — fijici e' -f Cj e'fjx = Xx — fijixj Ci Xi ef). 

Similarly, in the case of a link disconnection with dynamics 
reconfiguration we have 

^x = $a: -I- fij {ci e' -f ej e'- -I- e'- -f Cj e'f) x 

= Ax -I- {xi — Xj) Ci (xj — Xi) ej 


from which Theorem 2 follows. As for the case of a node 
disconnection without dynamics reconfiguration, we have 

$x = <I>x — fij (ci e'j 4- Cj e'i) x 

j&Mi 

— Ax I ^ ^ fij Xj I ei Xj ^ ^ fij Cj 

V'eAi / 

which leads to the statement of Theorem 3. Finally, if we 
consider the case of a node disconnection with dynamics 
reconfiguration, we have 

$x = <I>x + (cj e'i -\- ej e'j -f Cj e'j -\- Cj e'^j x 

j&Mi 

= Ax - ^ (xi - Xj) e, + ^ {xj - Xi) ej 

\j&Mi J j&Mi 

and, hence, condition Gil. Concerning the second part of 
Theorem 4, notice that when $ = — L then the i-th row of the 
condition $x = Ax can be written as 

l-A/j IXj 4“ ^ ^ Xj — A Xi . 

Thus, it can be seen that Xj = Xj for any j G A/i implies 
Axi = 0 which concludes the proof. ■ 

Proof of Theorem 5. Clearly, conditions (i) and (ii) are 
necessary because, otherwise, we would have non-null initial 
states leading to zero output trajectories and hence indis¬ 
cernible from the output. Further, condition (i) ensures that for 
any state x such that S"(A)x f 0 with A G spec($) \spec(T>) 
the term e^*MS{X)S'{X) related to the output evolution of 
{Q, $) is not null. Since this term is never present in the 
output evolution of {Q, 4>), we can conclude that {Q, <I>) is M- 
discernible from {Q, 4>) in such a state. Similarly, condition 
(ii) ensures that (0, $) is M-discernible from ((/,$) in all 
the states x such that S'{X)x 0 for at least one A G 
spec(4>) \ spec(4)). Hence, under conditions (i) and (ii), a 
necessary condition for a pair of states (x,x) to be M- 
indiscernible is that S"(A)x = 0 and S'{X)x = 0 for any A G 
spec(4))\spec(4>) and A G spec(4>)\spec(4)), or, equivalently, 
that X belongs to UAespec(<i.)nspec( 4 .)l^(A) and x belongs to 
UAespec(<i>)nspec(<i>)^(^)- Consider now two such states x and 
X. Condition (iii) ensures that, for any A G spec($)nspec(4>), 
one has MS{X)S'{X)x = MS{X)S'{X)x if and only if 
S{X)S'{X)x = S{X)S'{X)x. This, in turn, implies that {x,x) 
are M-indiscernible if and only if x = x G X, which proves 
the sufficiency of conditions (i)-(iii). Finally, the necessity of 
condition (iii) readily follows from the fact that, when such a 
condition is not satisfied for some A G spec(4)) fl spec(4>), 
there exist pairs of states (x, x) ^ Xp with x G V (A), 
X G V(A) such that MS{X)S'{X)x = MS{X)S'{X)x. m 
Proof of Lemma 1. Consider the joint observability matrix 
[On — On] which coincides with the standard sampled- 
data observability matrix associated with the pair (F, A). 
Since spec(A) = spec(<I>) U spec($) and A has dimension 
2n X 2n, we can apply the Kalman-Bertram criterion for the 


observability of sampled-data systems and conclude that, when 
conditions dnii and N > 2n are satisfied, the null space of 
[On — On] coincides with the set of unobservable states of 
the pair (F, A), which is precisely the set I{M). Hence, we 
have Onx = Onx if and only if (x,x) € X{M). Then, the 
statement of Lemma 1 can be proven by proceeding as in the 
proof of Theorem 2 of [25], to which the reader is referred 
for details. ■ 
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